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Critical temperature of the superfluid transition in bose liquids 
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A phenomenological criterion for the superfluid transition is proposed, wfiidi is similar to the Lin- 
demann criterion for the crystal melting. Then we derive a new formula for the critical temperature, 
relating Tx to the mean kinetic energy per particle above the transition. The suppression of the 
critical temperature in a sufficiently dense liquid is described as a result of the quantum decoherence 
phenomenon. The theory can account for the observed dependence of Tx on density in liquid helium 
and results in an estimate Ta ~ 1.1 K for molecular hydrogen. 



PACS numbers: 67.40.-w, 05.30.Jp, 67.40. Kh 

In connection with the search for new bose-condense 
systems, attention was again attracted recently to the 
possible superfluidity in molecular hydrogen H2 [y|] , dis- 
cussed earlier in |^-g. The very possibility crucially 
depends on the expected value of the critical tempera- 
ture Tx, since normally molecular hydrogen crystallizes 
at ~ 14 K and some measures should be taken to keep it 
in liquid phase at lower temperatures, which may appear 
harder to achieve if Tx is too low. Hence one needs first to 
evaluate T\ in supercooled hydrogen. In a more general 
context we are faced with the problem of determining Tx 
in a given quantum liquid. Usually, the first step is to use 
the well known formula for the Bose-Einstein condensa- 
tion temperature in the ideal gas 
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where n is the density and m is the particle mass {kb = 
1). Though this formula results e.g. in a rather reason- 
able estimate of Tx in liquid helium (^^ 3 K instead of 
the correct value 2.17 K at the saturated vapor pressure) 
it still seems unsatisfactory in case of dense liquids be- 
cause it cannot account for the observable dependence of 
Tx on density. According to (|]j) the critical temperature 
increases as n'^'^, while, on the contrary, Tx in liquid ^He 
slightly decreases when the system is compressed. Cer- 
tainly, one can argue that it is not the bare mass that 
enters in (|l|) but rather an effective one, depending on 
density Q, but there seems to be no simple and general 
expression for this effective mass. A calculation of Tx 
based on the Landau quasiparticle spectrum H can be 
applied only when this spectrum is already known. Qual- 
itatively the observed behavior of Tx can be explained 
also in the lattice model ||] but in general it is not so 
obvious how to relate the properties of a liquid to those 
of bosons on a lattice. For this reasons some new esti- 
mate of Tx is required, simple enough to serve as a first 
approximation, but which can account for nonmonotonic 
dependence of the critical temperature of a bose system 
on its density. In this Letter we propose such a new gen- 
eral formula, which relates Tx to the mean kinetic energy 
per particle in the normal phase — the quantity which 
behavior is now well understood [^-12|. 



The most general approach to superfluidity is, perhaps, 
to start with the Feynman path integral expression for 
the partition function Z for a system of A^ interacting 
bose particles M 



^=]^E/n^-'/n^r^M-p(-^^) 



(2) 
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where V(ri — Tj) is the interparticlc interaction potential 
and P = 1/T. The integration in (g) is over all paths 
with ri(0) = Ti, ri{hf3) = Pvi, where P is some permu- 
tation of A^ particles, and the sum in (0) is over all such 
permutations. 

At high temperatures only the identity permutation 
is important, since particles cannot move far away from 
their initial positions in 'time' hp. As /3 increases, a given 
path r(r) can spread on a larger distance, until suddenly 
it appears possible to end the path at the position of 
a neighboring particle. Then, as discussed in detail in 
the recent review 



y 



the rings of exchanges of arbi- 
trary length are formed, which can be shown to lead to 
the superfluid behavior. Let us now consider a system of 
distinguishable particles, take arbitrary particle and eval- 
uate its mean square displacement from initial position 
(((5r) ) in imaginary time. To estimate the critical tem- 
perature, when exchanges can no longer be neglected, we 
propose a criterion 
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where ^ is some numerical factor, to be determined later. 
This condition merely states that the mean displacement 
of a given particle in imaginary time is comparable to the 
interparticle spacing. The condition (0) is similar in spirit 
to the well known Lindemann criterion in case of crystal 
melting and is inspired by the visual representation of 
the paths, arising from numerical simulations [l3. 

In the ideal gas the left hand side of (0) is essen- 
tially the square of the de Broglie thermal wavelength 



A|, = h^ /mT, but interactions will tend to reduce {{5v) ) 
(this was observed e.g. in p4|). There are two main 
mechanisms of such a reduction. The first one is re- 
lated to decoherence due to interaction with environ- 
ment. Neighboring particles, in a sense, 'measure' the 
position of the particle we are looking at, thus reduc- 
ing its quantum uncertainty in coordinate space. Such 
a decoherence phenomenon was discussed a lot for more 
than two decades with respect to the transition from the 
quantum behavior to the classical one (see e.g. Ilq|). The 
second mechanism is more typical for crystals or glasses, 
where particles are almost localized by potential barriers. 
The problem of estimating of {{5v) ) in a system of 
interacting particles is still a very complicated one, even 
if exchanges are neglected. For a liquid not too close 
to crystallization it can be significantly simplified, how- 
ever, by treating the rest of the system as some sim- 
ple fluctuating environment. The simplest choice is the 
Caldeira-Leggett model |1^ , which describes interaction 
of a particle with a thermal bath of harmonic oscillators. 
In this model the particle motion in the imaginary time 
is governed by the effective action 
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where the kernel K{t) is determined by 
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where ujn — 2T:n/hp and I{uj) is the spectral density of 
bath oscillators (see e.g. [0). In our case one can view 
(Q) as a trial action and try to evaluate K{t) variation- 
ally, but it is easier to directly relate the parameters of 
the effective action to some observables of the system (see 
below). For quadratic action (m in three dimensions the 
mean displacement is 
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At small frequencies it seems natural to expect the 
ohmic behavior of the kernel K{t), corresponding to 
the linear friction, when damping is proportional to the 
velocity of the particle (cf |]l8[). This is modeled by 
I{uj) = ^muj, where 7 is the damping parameter, and 
Cn = 7|w„|. The Caldeira-Leggett model with such a 
dissipation kernel is a quantum analog of the standard 
Langevin equation with the friction 7717 and with the 
white noise random force. Then the sum in (0) is carried 
out and 
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where C ~ 0.577 ... is Euler's constant and ipix) is the 
psi function (the logarithmic derivative of the gamma 
function) . 

Now we have only one parameter, 7, which describes 
interactions in the system. The ideal gas limit is recov- 
ered at 7 ^- 0, when {{5y) ) ^ fi l3/2m. Then the cri- 
terion (|) resuhs in Tx = {l/2i)h^n'^/^ /m. Since this is 
the ideal gas formula for the critical temperature (H) we 
conclude that 1/2^ ~ 3.31, i.e. 



^-0.15 



(8) 



Though the criterion yj^ is useful only if ^ is some univer- 
sal constant we cannot exclude some weak dependence of 
^ on density. We expect that ^ may be slightly smaller 
than (g|) in a system with short range order since it is 
more likely for a particle in a liquid to find a neighbor 
for exchange at a suitable distance than in the ideal gas, 
where density fluctuations arc more important [|l9| . 

In the opposite case, when ?i7/3/27r ^ 1 the mean 
displacement diverges logarithmically i.e. ((<5r) ) ^ 
ln(?i7/3/27r). Though {{5v) ) still tends to infinity as 
T -^ 0, it is much less, than the de Broglie thermal 
wavelength of a free particle. This logarithmic behav- 
ior is related to the ohmic spectrum at small frequencies. 
Then for the critical temperature at 7 — > 00 we have a 
very simple formula 



Ta = a— exp(-C- — ) , To = —77/-' 
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(9) 



where a = exp(C). 

We see now, that the temperature of the superfiuid 
transition crucially depends on the ratio hrj JTq. At 
?i7 ^ To the estimate (||) is valid and the critical temper- 
ature is essentially Tq ^ r?^'^ , while at fi'^ ^ Tq the criti- 
cal temperature of the A transition is exponentially small 
due to the decoherence phenomenon. Qualitatively this 
can be understood as follows: in the process of exchange 
particles move through the viscous media and loose co- 
herence, needed for superfluidity to establish. Since 7 
should increase with increasing density, the formula ob- 
tained do describe the suppression of T\ in sufficiently 
dense systems. 

The approximation of a constant friction would be cor- 
rect if we dealt with a heavy Brownian particle, which 
moves slowly than particles in a liquid. In real liquids 
the ohmic spectrum /(w) ~ a; is physical only at fre- 
quencies lower than the collision rate and must have a 
cutoff at some frequency Wc ~ 7, which now should be 
taken into account. One can take e.g. the Drude model 
for the damping io 



C„ = 7a;^|w„|/(wc -I- |w„|) 



(10) 



Apart from the dispersion of the friction coefficient 
(memory effects) this imphes that the random force act- 
ing on a particle is correlated for times less than 1/wc- 

Qualitatively, however, the picture outlined above re- 
mains unchanged. If we e.g. put Wc = 7 then again 
the formula (O) is recovered but now with a = cxp(C + 
V37r/9). 

At low temperatures 7 can be related to the mean 
kinetic energy per particle. The kinetic energy in the 
Caldeira-Leggett model may be written as 



K{T)^-T- 
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The sum is easily evaluated for a model (ll^) and e.g. at 
Wc — ^ cxD (purely ohmic limit) we obtain 
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This very expression was derived earlier p(| for liquid 
■*Hc above the A point by taking the velocity autocorre- 
lation function to be of a simple exponential form with 
damping 7 (denoted in llfl] by o^o), just as for the Brow- 
nian particle. At T ^ ?i7 the kinetic energy tends to 
|r, while at low temperatures K{T) has a finite limit, 
denoted hereafter by K. This is quite a general behavior, 
valid in the Debye model as well ||21[] . Strictly speaking, 
K is the kinetic energy just above T\, since exchanges 
are not taken into account, but e.g. in helium the differ- 
ence between ii' ~ 16 K and the real zero-point energy 
is ~ 1.5 K and will be neglected here. For the particular 
model ( [I4 ) we have K = {3h'y/2n) InUjc/j), i-e., as was 
already mentioned in |10| (see also [gO|), the frequency 
7 up to a logarithmic factor coincides with K. Thus we 
conclude that in general case at low temperatures 



hj r^ K 



(13) 



For different I{lo) the proportionality coefficient is actu- 
ally cutoff dependent and e.g. for the model ( p^ with 
Wc = 7 one has Tij = \/3K. Since in the formula (0) 
the factor a also depends on a high frequency behavior 
of I{lu) we can finally write 



T^=AKexpi-B^) 



(14) 



for K ^ Tq, where A and B are some model dependent 
constants. For the Caldeira-Leggett model with damping 
kernel (|l^) at Wj, = 7 we have 

A = (\/3/27r) exp(C -I- \/37r/9) ~ 0.899, 

B = (V37r/6)C~ 0.907 C (15) 

The physical meaning of (H) is clear: if K, which at 
T <^ K may be viewed as an effective 'internal' temper- 
ature of the system pffl, is much larger than the tran- 
sition temperature in the ideal gas, the superfluidity is 
suppressed. 



Let us now compare the formula (14) with experimen- 
tal data for liquid helium. We need then some explicit ex- 
pression for K{n). There are different estimates of zero- 
point energy in helium (see e.g. |22|JlC| ]). All of them are 
in general consistent both with the experimental data pl| ] 
and with results of Path Integral Monte Carlo (PIMC) 
calculations P,f2|. The kinetic energy increases with 
density due to the repulsion core in interatomic potential. 
Here we shall use the London's formula 



K = 



2'Kn^d 



(16) 



m(a-0.891d)2(a-H 0.713d) 
where a — n"^/^, d ~ 2.4 A |22|. This formula is sim 



pie and transparent being an interpolation between low 
density limit K ~ h^dn/m of the energy of a gas of hard 
spheres of radius d, and a quantum mechanical estimate 
K ^ h /mia — dof' with dp ~ d in the high density limit. 




Density ( nr^,'^ ) 

FIG. 1. Temperature of the superfluid transition vs the 
reduced density n* — nr^ (ro = 2.556 A). Solid line is the 
theory at ^ = 0.12, crosses denotes experimental data for 
helium, dashed line corresponds to the ideal gas. 

With this expression for K the formula (^^ with co- 
efficients (|l5|) at C = 0.15 results in Ta ~ 1.35 K, which 
is lower than the experimental value of 2.17 K. We may 
recall, however, that ^ may be smaller in dense system 
than in the ideal gas, and try to fit ( [l4| ) to experimental 
data treating ^ as an adjustible parameter. The result 
of such a fit is shown in Fig. 1. Here the critical tem- 
perature is shown as a function of the reduced density 
n* = n^Q, where tq — 2.556 A is the length parameter 
of the Lcnnard- Jones interatomic potential for helium. 
Experimental values of T\ for liquid helium, taken from 
the very accurate empirical expression for the A- line p3[ , 
are shown by crosses. The dashed curve represents the 
ideal gas formula (Q) , while the solid one is the best fit of 
T\ from (n^ and ( |15D to experiment, which corresponds 
to ^ ~ 0.12. This value is only slightly smaller than the 
ideal gas limit ^ = 0.15. Given the simplicity of the as- 
sumptions, the agreement of the theory with experiment 



is quite satisfactory. Since BK/Tq ~ 1 at n* ~ 0.2, we 
can not expect equation ( p^ to be quantitatively valid 
at such a low density, but still it is clear that the theory 
really can describe the crossover from the ideal gas be- 
havior to the observed dependence of Tx on the density. 

Now, as far as the molecular hydrogen is concerned, 
we may take the London's formula (Iq) with d ~ 2.7 A 
p4| , which corresponds to a stronger interparticle poten- 
tial. Assuming the numerical coefficients A and B to be 
the same as in helium (with ^ ~ 0.12) we obtain from 
( p^ ) Ta ~ 1.1 K for density n ~ 26 nm"^. Even if we 
reduced the density to that of helium i.e. n ~ 22 nm^^ 
the critical temperature would be only ~ 2.1 K. This is 
much less optimistic than the original estimate 6 -^ 8 K 
Q , based on the ideal gas formula (|l|) , though consistent 
with later estimates [Q and with the PIMC analysis of 
finite hydrogen clusters [g5[. The value ~ 1 K, obtained 
recently for hydrogen films with impurities ||l| is also of 
the same order. 

In conclusion, we have obtained a general estimate 
for the critical temperature of the superfluid transition 
in a bose liquid. Starting from the phenomenological 
Lindemann-like criterion (0) for the transition and mod- 
elling the decoherence effect, which suppresses T\ in liq- 
uid, by the Caldeira-Leggett model, we arrive at a sim- 
ple expression (^, relating T\ to the kinetic energy K 
(which is essentially the zero-point energy per particle). 
The fast increase of K with density due to the repulsion 
core in interatomic potential accounts then for the reduc- 
tion of Tx in dense systems, when K '3> Tq ^ h n^/'^/m, 
as observed e.g. in liquid helium. Using the London's in- 
terpolation formula for K{n) and one fitting parameter 
the model can even quantitatively describe experimental 
behavior of T\ . 

Our approach suggests that, in a liquid, contrary to 
weakly nonideal gas, Tq may not correspond to any char- 
acteristic temperature. The temperature where quantum 
effects become important is of the order of -ftT, and is much 
larger than Tq [0, while T\ at high densities is smaller 
than To, and the difference increases with density. The 
error of using the ideal gas formula (|l|) is not so much for 
helium, but may be of importance in molecular hydrogen 
(which is more dense, due to stronger interaction) , where 
our formula gives much smaller value of T\ . 

An important question, which remains open within the 
phenomenological approach, concerns the universality of 
the numerical constants in (nj ). Though the parameter 
^ (analogous to the Lindcmann ratio in case of crystal 
melting) is to some extent fixed by the ideal gas limit (0) 
further work is needed to clarify its possible dependence 
on system parameters. 
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